Geometric approach to asymptotic expansion of Feynman integrals 
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We present an algorithm that reveals relevant contributions in non-threshold-type asymptotic 
expansion of Feynman integrals about a small parameter. It is shown that the problem reduces to 
finding a convex hull of a set of points in a multidimensional vector space. 



I. INTRODUCTION 

O ■ „ , 

' Evaluation of Feynman integrals depending on multiple parameters is a notoriously difficult task. When direct 
computation fails, one resorts to studying asymptotics in various limits. In practice, a few first terms in the expansion 
may already suffice to reach the desired precision. However, expansion of a multi-loop integral may become non- 
trivial due to an interplay of parameters with the integration variables (components of loop momenta) . Classification 
Q \ of relevant sectors in the integration space is itself a challenging problem P, [3] • 

One important case is the asymptotic expansion in momenta and masses in the limits typical for the Euclidean 
space. This problem has been completely solved in terms of sums over subgraphs 0-0] • At least one automated 
tool [H, Q implements this approach in practice. For a more general situation, including the limits appearin g in the 
Minkowski space, there exists the universal strategy of expansion by regions in the momentum space fll. 13, lldf. In 
1— 1, all known cases, it produces correct results, but a rigorous proof is still lacking. Typically, one manually analyzes a 
r& ' multi-scale problem, starting from simpler examples that can be checked against known analytical results or numerical 



Oh 1 estimates, computed e.g. with FIESTA (later versions [l2| of FIESTA may also evaluate a few first terms in a 
fS ■ given asymptotic expansion.) 

An important type of non- Euclidean expansions, the so-called threshold expansion requires the most careful 
treatment. Cancellation of dominant terms becomes obvious only in a specially chosen frame or with a certain routing 
of loop momenta. In what follows we try to elaborate some approach to non-threshold asymptotic expansion, based 
, on alpha-representation of integrals, and describe a simple practical algorithm. 
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\Q , II. EXPANSION BY REGIONS AND ALPHA-REPRESENTATION 

00 

A thorough introduction to the expansion by regions and alpha-representation can be found elsewhere [Tol [l3| . 
Here we briefly introduce the basic notation with a trivial example. 

Consider a family of one-loop propagator-type integrals in the Euclidean space: 

hia^az-p 2 ,™?) = J AA__J_ D 1= k 2 +m 2 , D 2 = (k + p) 2 + m 2 . (1) 

A specific integral is determined by the exponents a\ and a 2 and depends on the two parameters, m 2 and p 2 . The 
structure of the expansion does not depend on a\ and a 2 and we will not mention those exponents in the following 
" ' discussion. 

We consider the asymptotics of Ii(ai,a 2 ;p 2 ,m 2 ) in the limit when \p 2 \ ^ m 2 , or p = \m 2 /p 2 \ <C 1. The naive 
Taylor expansion does not capture the complete asymptotic behaviour since the integration variables (components of 
k) span all values from — 00 to +00, and in particular can be as small as to or as large as y[p^|. 

The prescription in this case is to find regions, or scalings of momentum components that after the expansion 
provide non-zero contributions. In each region, we first Taylor expand the integrand and drop the scaling restrictions. 
In our example, there are three non-zero regions, summarized in Tab. [I] For example, in the region (c) one expands 
D 2 as follows: 



h{ai,a 2 ;p 2 ,m 2 ) 



d D k 1 
(2ir) D (k 2 + m 2 ) a i 



1 a 2 (k 2 + 2kp + to 2 ) 



(p2)a 2 (j,2)a 2 + l 

d D k 1 f d D k (...) 



{2tt) d (k 2 +TO 2 ) Q i(p 2 ) a 2 J (2tt) d (k 2 + m 2 )^(p 2 )^ +1 



other regions (2) 
- other regions. 



In the last line, we dropped the restriction on k and the problem reduced to evaluation of (multiple) Feynman 
integrals with simpler denominator factors, various denominator exponents and possibly more complex numerators. 
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Di = k 2 + m 2 , 
£> 2 = (fc+p) 2 +m 2 , 
p = \m 2 /p 2 \ < 1 


W = Si + IE2, 

J 7 = a;iX2 (p 2 + 2m 2 ) + x 2 m 2 + xlm 2 


.' k ' 

p : 
> ■ 


(b) 


fc 2 ~ \p 2 \ > m 2 , 
D< 6) = fc 2 , Di b) = (k + p) 2 


a?i, X2 ~ p A , 

7 /(&) T-(6) 

U — XI + X2, -T — XlX2 




(c) 


\k 2 \ ~m 2 , 
Z?< c) = k 2 +m 2 , D^=p 2 


A A 1 

£2 ~ P , Xl ~ p , 

7 /(c) t-(c) 2 . 2 2 
W = Xl, J" = X1X2P + X]Tn 





(d) 


|(fc+p) 2 | ~m 2 , 
D[ d) = p 2 , £>W = (fc + p) 2 +m 2 


A 1 A 

a;2 ~ p , £i ~ p 

7 /(d) T-(d) 2 2 2 
W v ' — X2, J 7 — X1X2P + x 2 m 



TABLE I: Regions (b-d) of expansion of a double-scale integral (a) in the momentum space and in the alpha-representation. 



The non-trivial statement is that the double-counting which could have been introduced disappears in the sum of 
all regions. For the purpose of the following discussion we assume that the tensor reduction of numerators and the 
proliferation of terms can be managed; we will focus on the transformations of denominator factors in every region 
(e.g. (k + p) 2 + m 2 — > p 2 in the example above). 

Some types of asymptotic expansion may require more elaborate choice of regions. For example, fco may scale 
differently from fcj or some combination of components may have a separate scale (this happens, e.g., in Sudakov 
limits). It is thus desirable to have an explicitly covariant formalism to identify regions, independent of the frame 
choice and the routing of momenta. For that purpose, we may switch to the alpha-representation of integrals. We 
re- write an integral with n lines (denominator factors) over Z?-dimensional loop momenta as an integral over n positive 
parameters x\, ..,x n . Information about the graph is then encoded in the two homogeneous polynomials, IA and J- . 
For example, our integral above is 

h{a u a 2 -p 2 ,m 2 ) = ^ ± °* |~ D/2) [°° d Xl dx 2 5(l - x x - x 2 )x^~ x x^~ x W^- D T° , (3) 

r{a 1 )T(a 2 ) Jo 

U = x\ + X2, T = xiX2(p 2 + 2m 2 ) + x\m 2 + x 2 m 2 . 

The expansion by regions may also be formulated in the alpha-representation 0, [Tol | . Instead of finding the scaling 
behaviour of loop momentum components, here we deal with the scaling of each parameter x% that directly corresponds 
to the scale of the i-th line (denominator factor) of the original integral. During the expansion, only the leading terms 
remain in the polynomials U and J 7 , and the resulting alpha-representation represents the integrals obtained by the 
expansion in the momentum space. The last column of Tab. |T] demonstrates the scaling of alpha-parameters and the 
polynomials corresponding to each region. 

Note that in the language of alpha-parameters, the difference between the threshold-type and non-threshold ex- 
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pansion becomes clear. Let us consider integral of Eq. Q]in the threshold limit y = m 2 + ^- <C m (that, of course, 
implies that p 2 < 0, i.e. this limit is essentially non-Euclidean). Choosing the frame where p — (po,0) and re-routing 
the loop momentum, we obtain the denominator factors D\ = k 2 , + k 2 + kgpg + y and D 2 = k 2 , + k 2 — fcoPo + y- 

This integral has two non- vanishing regions. The first "hard" region is characterized by k <~ m. The second 
"potential" region corresponds to kg ~ y/m, \k\ ~ In the language of alpha-parameters, in the "hard" region 

only the second term survives in the polynomial T — y{x\ + X2) 2 — ^r(xi — x 2 ) 2 - The most troublesome "potential" 
region stems from a thin layer in the integration space near the surface x\ = x 2 , when the second term has the same 
scaling as the first. 

In a similar way, more complex threshold expansions receive contributions which depend on cancellations between 
the terms in the expanded J- which happens along some non-trivial surface and not at zero or infinity. Presently we 
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do not know a general rule to identify such surfaces and find substitutions revealing such regions. Instead, we focus 
on the "usual" regions that can be determined by examining independently the monomials in U and J- ' . However 
limited, this problem is still important for many applications. 



III. GENERAL FORMALISM 

We consider an Z-loop Feynman integral 

d D k 1 ...d D ki 

{2n) lD D a 1 1 ...D^ 1 ( ' 

which depends on n exponents a\,...,a n , scalar products of e external momenta p\,...,p e and parameters (such as 
masses) in denominator factors Di. The latter must be quadratic in momenta but other than that may have any 
form, e.g. correspond to a massive, such as —{k% + Pj) 2 + m 2 — iO, or a static propagator, such as (—2kiPj ± iO). The 
alpha-representation for this integral has a general structure 

I(ai,...,a n ) = c dx\...dx n 8{1 — x\ — ... — x^x^ 1 ^ 1 ...x'^ l n ~ 1 U a J 7b , (5) 
Jo 

where coefficient c and exponents a and b depend only on I, D, and at. U and T are homogeneous polynomials (of 
order I and I + 1, respectively) of integration variables Xi, and J- also depends on the kinematic invariants. If the 
denominators Di correspond to some graph and have a standard form — k 2 + m 2 — iO (in Minkowski space), then the 
functions U and J- can be read off the graph in terms of trees and 2-trees 13]. In a more general case, one may obtain 
IA and J- with a tool found at http://www-ttp.particle.uni-karlsruhe.de/~asmirnov/Tools-UF.htm. In what 
follows, we will only discuss the properties of U and T that are independent of specific indices ctj. 

In dimensional regularization, "scaleless" integrals (having no inherent scale) turn to zero. More specifically, an 
integral is scaleless if it is possible to re-scale some loop momenta or their components so that the result remains 
proportional to the original integral, or Di({kj}, {aki}) — a Ui £),({&}), with some subset {ki} of integration momenta. 
In particular, massless vacuum bubbles vanish: 

/= f^= t dD ^ =a°-*>I = (6) 

In terms of the alpha-representation Eq. [5j a similar statement applies to homogeneity of U and T with respect 
to a subset {B} of integration variables Xi ({B} should not coincide with the full set of {xi}). Integrals vanish if 
U({xj}, {axi}) — a u U({x\) and F{{xj}, {axi}) = a' ^({x}), i £ {B}, with some scaling dimensions u and /. 

In order to avoid separate treatment of U and J 7 , one may consider the product UT that incorporates the scaling 
and asymptotic properties of both factors (but may contain many terms) . 




IV. GEOMETRIC INTERPRETATION OF ASYMPTOTIC EXPANSION 



Let us start with some integral in the alpha-representation Eq. [5] with integration variables xi,...,x n and a small 
expansion parameter p. Each of M terms in T corresponds to a vector of n + 1 exponents (we here neglect common 
factors and numeric coefficients, irrelevant to the non-threshold expansion): 

p r °x T ^ ...x r ™ -> (r ,ri,...,r„), (7) 

and J- corresponds to a set {F} of M points in (n + l)-dimensional vector space. Due to homogeneity of J 7 , all these 
points belong to an n-dimensional hyperplane n + ... + r n = I + 1, parallel to the 0-th axis (the axis of tq). 

Terms of U have no explicit powers of p in the coefficients. The corresponding set {U} is thus confined to an 
(n— 2)-dimensional hyperplane ro = 0, r% + ... + r n = I. In Fig.[T]we present such points corresponding to the example 
in Eq. [TJ where the three terms of {F} are denoted with crossed points and the two terms of {U} with diamonds. 

If we fix the scales of alpha-parameters as Xi ~ p Vi , then the scale of a monomial is p r °x r l 1 ...a^/* ~ p r °+ v ii"i+...+r n u n ^ 
p rv with r — (tq, r n ) from {F} and v = (1, v\, v n ). Graphically, rv represents the length of a projection of the 
vector r on the direction v. 

Some special choices of directions v determine the regions of expansion that we seek. The terms in T that remain 
after the expansion are all characterized by the same scale in powers of p. All points of the corresponding subset {F'} 
then feature the same value of the projection on v, i.e. these points belong to the hyperplane orthogonal to v. 
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FIG. 1: Graphical representation of sets {F} (crossed points) and {U} (diamonds) corresponding to the integral of Eq. [T] 



The points corresponding to the neglected terms will be located "above" this hyperplane (since v always points 
"up" with respect to the 0-th axis). In other words, {F'} belong to a facet of the envelope, or the "convex hull" of 
the set {-F 1 }, while the corresponding v is the normal vector to that facet. In a similar manner we may define a subset 
{[/'} of the remaining terms in li. 

Relating the three expansion regions of Tab. [I] to graphics in Fig. [TJ we find the corresponding points and vectors 
(points as denoted in the figure): 

. (b): v= (1,0,0), {F'} = (C), {U>} = (D, E), 

. (c): v= (1,1,-1), {F'} = (A,C), {{/'} = (D), 

. (d): v = (1, -1, 1), {F'} = (B, C), {£/'} = (E). 

Here we exploit the freedom to re-scale all Xi by the same power of p, i.e. shift v by any vector a = (0, A, A). If 
v corresponds to a region, then v' — v + a determines the same region. For example, v' = (1, 2, 0) also corresponds to 
the region (c) above. It is convenient to choose v parallel to the plane where points {F} are confined, i.e. orthogonal 
to the vector (0, 1, 1). 

Normally, only a few scaling choices produce non-zero regions. In our example, the choice X\ ~ p 2 , x-i ~ p°, leading 
to U = x 2l F = x^m 2 , or {F'} — (A), {U'} = (D), corresponds to a scaleless integral. As discussed above, this 
implies an existence of a scaling leaving both U and F invariant up to a pre- factor (in this case, x\ — » ax\). 

The requirement that a region does not vanish can be easily formulated in the geometrical language. Consider the 
polynomial UF and the corresponding set of points {UF}. After the expansion with the chosen scalings, we are left 
with its subset {UF'}. Numeric coefficients and kinematic invariants are irrelevant to the scalefulness of the region, 
and we get rid of them by projecting {UF'} on the plane r = 0. The thus obtained set of points {UFq} belongs to 
the (n — l)-dimensional hyperplane. 

Scalelessness implies that all terms of the polynomial UF' are homogeneous with respect to a certain re-scaling. We 
thus deduce that the points of {UFq} must then belong to an orthogonal space of the corresponding re-scaling vector 
iJh. In other words, {U Fq} is confined to at most (n — 2)-dimensional subspace, and its (n — l)-dimensional volume is 
zero. The latter property can be easily checked (and used to check whether a given integral vanishes). However, it is 
easy to see that the "bottom" facets of the convex hull for {UF} automatically correspond to scaleful regions: their 
dimension is (n — 1) by construction (otherwise they become "ridges" or "vertices"), and they (by selection) are not 
orthogonal to the plane ro = (thus the projection has non-zero volume). 

Finally, we may formulate the general procedure to determine the expansion regions. We start by building the 
set of points {UF}. Next, we find the n-dimensional convex hull C of the set {UF} in the n-dimensional plane 
ri + ... + r n = I + 1, using any preferred algorithm. The implementation that we chose, QHull does not allow 
building hulls of dimensionality lower than the dimension of vector space. Thus, one has to introduce local coordinate 
system and deal with non-integer coordinate values. However, it is also possible to project {UF} along any of axes 
ri, i 0, e.g. consider (n — l)-dimensional points rji = (ro, rx, r n _i). Convex hull C built for this projection will 
be the projection of the "true" convex hull C. Its dimensions will be stretched but the correspondence of the points 
to the facets and the vertices will persist. 
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FIG. 2: Double-scale two-loop vertex integral. 

From the (n — l)-dimensional facets of C we then select the "bottom", i.e. facets with normal vectors v pointing 
"up" , with v > 0. For each of those "bottom" facets, we choose the normal vector v such that vq = 1. Its components 
1 to n represent the relative scales of alpha-parameters x% and thus uniquely determine an expansion region. 

V. LESS TRIVIAL EXAMPLE 

Let us consider the integral in Fig. [5J this time defined in the Minkowski space (this example was first considered 
in 0] and [13, Chapter 10): 

2 2^ f d D k 1 d D k 2 
h( ai ,...,a 6 ;s,m ,M ) = J ^yw jy*y _jy* > ( 8 ) 

D x = (pi-fci-fc 2 ) 2 -M 2 , D 2 = ( Pl -k 2 ) 2 -M 2 , D 3 = (p 2 + fc a + fc 2 ) 2 - to 2 , 

Di = (p 2 + k 2 ) 2 -m 2 , D 5 = k 2 , D 6 = k 2 , 

p\ = M 2 , p\ = rn 2 , (pi + p 2 ) 2 = s, s, M 2 > m 2 . 

With S = to 2 + M 2 — s = —2pip 2 , its alpha-representation reads: 

U = x\x% + X3X2 + X5X2 + xiXi + X3X4 + X!X 5 + £3X5 -I- X4X 5 + xiXq + x 3 x 6 + x 5 x e , (9) 

T = M 2 x\x 2 + M 2 x\x A + M 2 xjx 5 + M 2 xix\ + M 2 x\x-s + M 2 x\x b + M 2 x\x & + 2M 2 xix 2 x 5 

+ m 2 x 2 x\ + m 2 x 3 X4 + m 2 X\x\ + rr?x^x\ -I- to 2 x 2 2;5 + m 2 x\x^ + m 2 x 2 xe + 2m 2 X3X4X5 

+ SxiX 2 X3 + SxiX 2 X4 + SX1X3X4 + SX1X3X5 + SX1X4X5 + SX1X3XQ + Sx 2 Xt,X4 + Sx 2 X3X$ + Sx 2 X4X§. 

For simplicity, let us analyze T instead of the product UJ-. Choosing to as the small parameter and preserving the 
order of terms, in the 7-dimensional space we have 25 points: {F} = (0,2,1,0,0,0,0), (0,2,0,0,1,0,0), (0,2,0,0,0,1,0), 

The projection of {F} along the 6-th axis (i.e. {F} p — (0,2,1,0,0,0), (0,2,0,0,1,0), ...) has a six-dimensional 
convex hull with 18 facets. Of them, four belong to the "bottom". Restoring the 7-dimensional normal vectors with 
unit 0-th component, we find: v x = (1,0,-2,-2,-2,-2,-4), v 2 = (1,0,0,-2,-2,-2,-2), v 3 = (1,0,0,-2,0,0,0), 
V4 = (1,0,0,0,-2,0,-2). Since U is scaleful, we also have to add the "hard" region: vq = (1,0,0,0,0,0,0). (The 
latter would appear automatically, had we used the product UF.) 

For illustration, let us consider the most non-trivial "ultrasoft-collinear" region corresponding to v\, or the scaling 
of alpha-parameters x\ ~ to , x 2 ~ 1/to 2 , X3 ~ 1/to 2 , X4 ~ 1/to 2 , x§ ~ 1/to 2 , and xq ~ 1/m 4 . First, we may check 
that vi is indeed orthogonal to the plane containing the points 5, 6, 15, 22, 23, 24, and 25 from {F} (in the order as 
in Eq. [SJ. Those points correspond to the terms remaining in F' after the expansion: M 2 x\x3 + M 2 x\x^ + . . .. 

In the momentum space, the interpretation becomes clear only in the special reference frame, where p\ = (M, 0, 0), 
p 2 = Mn + + Ti—, and n± — (1/2, 0, =pl/2). We also decompose the first loop momentum in plus- and minus- 

parts, fci = (fc + + k-,k, k + — To reproduce the "ultrasoft-collinear" region, we should prescribe the following 
scales to the components of loop momenta: k + ~ to 2 /M, fc_ ~ M, k ~ to, k 2 ~ m 2 /M. 

After the expansion, the denominator factors scale as: Dj ~ M 2 , D 2 , D3, D4, D$ ~ to 2 , D 6 ~ m A /M 2 . One can 
easily see how the powers of to here correspond to the components of — v±. 

VI. IMPLEMENTATION 

We wrote a Mathematica program that determines the expansion regions of a given Feynman integral based on the 
procedure described above. The general problem of building a convex hull of M points in n dimensions is well-known 
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in computational geometry; we employ the algorithm quickhull [l4[ that has complexity O (ML d / 2 J). It is sufficient 
when the number of lines is not too large; for example, finding 11 expansion regions of a 4- loop integral with 10 lines 
takes about 10 seconds on a laptop PC. 

The program has been checked against some non-trivial examples discussed in (loj and [l5j |. The code can be 
downloaded from http : //www-ttp . particle . uni-karlsruhe . de/~asmirnov/Tools-Regions . htm 

In order to run the program, one has to install the open-source package QHull [14]. If the executable is not in the 
current directory, Options [QHull] must be updated in the file asy.m. The program is loaded with command <<asy .m. 
The main function is AlphaRepExpand[ks,ds,cs,hi] , where ks is the list of loop momenta (e.g., {vl,v2}), ds are 
the denominators (e.g., {(pl-vl-v2) ~2-M~2, (pl-v2) "2-M"2, (p2+vl+v2) ~2-m~2, (p2+v2) ~2-m~2,vl~2,v2~2}), cs 
contains the kinematic constraints (e.g., {pl~2->M~2 ,p2~2->m~2 ,pl*p2->-S/2}), and hi represents the scalings of 
kinematic invariants with respect to the small parameter x (e.g., {M->x~0,S->x~0,m->x~l}). 

The output is a list of vectors specifying the scales of the alpha-parameters factors, 
or the non-zero components of vectors V{. For the last example above, the output is 

{{0 ,-2,-2,-2,-2, -4} ,{0,0,-2,-2,-2, -2} , {0 , , -2 , , , 0} , {0 , , , -2 , , 2} , {0 , , , , , 0}}, correspond- 
ing to the regions (us-lc), (lc-h), (h-lc), (lc-lc), and (h-h). These regions can be understood by analogy to the 
example of Section [V] 

VII. CONCLUSION 

We present an algorithm to find the relevant regions of expansion for a Feynman integral in a given limit of momenta 
and masses. The algorithm is implemented in Wolfram Mathematica language and uses open-source package QHull. 
The program and examples can be downloaded from our web-page. In the future, we plan to extend the code in order 
to apply it to more general parametric integrals. 
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